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AUTOMATED  TORSION  PENDULUM: 
CONTROL  AND  DATA  COLLECTION /REDUCTION 
USING  A  DESKTOP  COMPUTER 

John  B.  Enns  and  John  K.  Gillham 
Polymer  Materials  Program 
Department  of  Chemical  Engineering 
Princeton  University 
Princeton,  New  Jersey  08544 


A  torsion  pendulum  Interfaced  with  a  desktop  computer 
form  an  automated  Instrument  for  dynamic  mechanical 
characterisation  of  polymeric  materials.  The  com¬ 
puter  controls  the  initiation  of  the  oscillations, 
collects  the  digitised  data  and  calculates  the  shear 
modulus  and  loss  modulus  from  the  dsmped  oscilla¬ 
tions,  utilising  one  of  four  methods  of  analysis: 

1)  fitting  the  data  points  about  the  maxima  and 
minima  to  a  quadratic  equation  to  obtain  their  times 
and  amplitudes,  from  which  the  frequency  and  log¬ 
arithmic  decrement  can  be  calculated;  2)  fitting 
the  data  to  a  four-parameter  equation  of  srotion  by  a 
least  squares  technique;  3)  fitting  the  data  to  a 
six-parameter  solution  to  the  equation  of  motion  by 
a  non-linear  least  squares  technique;  and  4)  taking 
the  Fourier  transform  of  the  data,  which  results  in 
a  maximum  at  the  frequency  of  the  oscillation  whose 
amplitude  is  inversely  proportional  to  the  damping 
coefficient.  The  advantages  and  disadvantages  of 
each  method  are  discussed  and  the  results  of  tor¬ 
sion  pendulum  and  torsion  braid  analysis  (TBA)  ex¬ 
periments  sre  compared. 


The  torsion  pendulum  has  proven  to  be  an  important  and  versa¬ 
tile  tool  in  the  study  of  dynamic  mechanical  properties  of  mater¬ 
ials.  In  our  laboratory  it  has  been  applied  primarily  to  poly¬ 
mers,  although  elsewhere  it  has  been  used  with  a  wide  variety  of 
materials,  ranging  from  liquids  to  metals  and  ceramics.  The  basis 
of  its  wide  appeal  lies  in  its  fundamental  simplicity:  informa¬ 
tion  about  the  complex  modulus  of  the  material  under  Investigation 
is  obtained  by  simply  observing  the  decaying  oscillations  of  the 
pendulum.  After  the  pendulum  Is  set  in  motion,  it  is  permitted  to 
oscillate  freely  at  its  resonant  frequency  while  the  amplitude  of 
the  oscillatory  wave  decays.  .In  an  unautomated  system  it  is  a 
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relatively  simple  but  tedious  task  to  calculate  the  shear  modulus 
and  the  loss  modulus  from  the  period  of  the  oscillation,  its  log¬ 
arithmic  decrement  and  the  geometric  constants  of  the  system.  The 
independent  variable  in  the  investigation  of  dynamic  mechanical 
properties  of  a  material  la  often  temperature,  but  it  can  also  be 
time,  as  in  the  case  of  chemically  Teactlve  or  physically  aging 
systems. 

A  variation  of  the  torsion  pendulum,  torsional  braid  analysis 
(TBA) ,  utilizes  a  supported  specimen  so  that  the  dynamic  mechani¬ 
cal  properties  of  a  sample  can  be  monitored  in  the  liquid  as  well 
as  the  solid  states  (.1,  2) .  An  inert  multlfllamented  glass  braid 
is  impregnated  with  the  Temple  (usually  in  its  liquid  state  or  in 
solution) .  The  observed  dynamic  mechanical  properties  are  rela¬ 
tive  due  to  the  composite  nature  and  complex  geometry  of  the  spec¬ 
imen. 

The  purpose  of  this  paper  is  to  describe  an  automated  torsion 
pendulum  controlled  by  a  desktop  computer,  to  discuss  four  separ¬ 
ate  methods  of  data  analysis,  and  to  compare  the  results  of  a  tor¬ 
sion  pendulum  experiment  and  a  TBA  experiment  using  the  same  epoxy 
resin. 

Instrumentation 


A  schematic  diagram  of  the  torsion  pendulum  is  shown  in  Fig¬ 
ure  1.  Free  oscillations  are  initiated  by  an  angular  step-dis¬ 
placement  of  the  upper  member  of  the  pendulum.  The  response  of 
the  lower  member  is  a  damped  wave  at  the  natural  frequency  of  the 
system,  and  therefore  is  related  to  the  physlco-mechanlcal  proper¬ 
ties  of  the  specimen. 

The  damped  oscillations  are  converted  to  an  electrical  signal 
by  a  non-drag  optical  transducer:  light  is  passed  through  a  pair 
of  polarizers,  one  of  which  serves  as  the  inertial  mass  of  the 
pendulum,  to  a  photo-detector.  The  temperature,  humidity  and  gas 
(usually  helium)  surrounding  the  specimen  are  closely  controlled. 

The  torsion  pendulum  has  been  interfaced  with  a  digital  desk¬ 
top  computer  (Hewlett  Packard  9825B)  shown  in  the  system  diagram 
Figure  2  (3).  The  motors  which  align  the  specimen  and  initiate 
the  waves  are  under  computer  control  via  the  scanner  (HP  3495A) 
and  relays.  At  present  the  direction  of  the  temperature  scan  and 
the  status  of  the  experiment  (whether  to  hold,  reverse,  or  termi¬ 
nate)  at  either  of  the  temperature  limits  set  by  the  programmer 
(Eurotherm  Corp.)  are  under  computer  control  as  well,  but  the  rate 
of  temperature  change  and  the  limits  are  not.  The  amplified 
thermocouple  and  wave  signals  are  digitized  by  a  high  speed  digi¬ 
tal  voltmeter  (HP  3437 A)  whose  scan  rate  is  programmable,  and  the 
scanner  supervises  the  X/0  activity.  The  computer  calculates  the 
frequency  and  damping  parameters  from  the  raw  data  and  plots  the 
dynamic  mechanical  properties  of  the  specimen  as  a  function  of 
temperature  and/or  time.  A  photograph  of  the  equipment  is  shown 
in  Figure  3.  A  commercial  version  of  the  automated  torsion 
pendulum/torsion  braid  analyser  is  available  from  Plastics  Analysis 
Instruments,  Inc.,  P.0.  Box  406,  Princeton,  Hew  Jersey. 


For  each  damped  wave  the  computer  goes  through  a  control 
sequence,  schematically  represented  in  Figure  4.  Since  the  spec¬ 
imen  may  twist  due  to  an  uneven  distribution  of  thermal  stresses, 
the  alignment  motor  rotates  the  pendulum  through  a  gear  train  to 
the  same  reference  position  at  the  start  of  each  control  sequence. 
To  initiate  the  oscillations,  a  second  motor  rotates  the  pendulum 
a  specified  angular  displacement  against  the  tension  of  a  spring. 
The  pendulum  is  held  in  this  cocked  position  until  oscillations 
set  up  by  the  alignment  and  cocking  procedure  have  decayed,  at 
which  time  the  clutch  is  disengaged  and  the  pendulum  swings  back 
so  as  to  oscillate  about  the  reference  position.  The  data  are 
then  collected  and  reduced.  The  temperature  (or  time,  for  iso¬ 
thermal  runs)  is  measured  with  the  specimen  in  the  cocked  position 
and  again  after  the  data  are  collected.  After  plotting  the  re¬ 
duced  data,  the  oscillation  is  monitored  until  it  decays  to  within 
specified  limits  and  the  cycle  repeats. 

The  data  obtained  from  the  torsion  pendulum  can  be  displayed 
in  various  modes  (4):  the  shear  modulus  G'  is  given  by 

G’  -  KI(-^M1  +  (1) 

or  by  its  approximation 

G»  a  4it2KI(^)2  (2) 

where  F  is  the  period,  4  is  the  logarithmic  decrement  and  K  is  a 
geometric  constant.  In  a  TBA  experiment,  where  K  is  unknown,  the 
relative  rigidity  t*(l/P)2]  is  measured.  Usually  the  logarithmic 
decrement  term  in  equation  (1)  is  negligible;  only  in  the  transi¬ 
tion  regions,  where  4  >  0.6,  does  it  become  greater  than  one  per¬ 
cent.  In  Figure  5  both  the  shear  modulus  and  its  approximation 
are  plotted  (5,  6):  the  curves  are  indistinguishable  except  in 
the  transition  regions.  The  energy  lost  during  the  deformation 
can  be  displayed  in  a  variety  of  ways  (Figure  5):  loss  modulus 

G"  -  4*KI»/P,  (3) 

logarithmic  decrement 

r»i 

4  «  »  jr  »  aP  ■  »tan$,  (4) 

and  damping  coefficient  a.  A  shift  Is  observed  in  the  maximum  by 
which  the  transition  temperature  is  identified:  for  example, 
Tg(G")  <  Tg(4)  <  Tg(o)  for  a  solid-to- rubber  transition.  (For  a 
rubber-to-solid  transition  the  shift  occurs  in  the  reverse  order.) 
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Software 


An  efficient  algorithm  is  required  to  monitor  the  oscillatory 
wave  signal  in  real  time.  The  flow  chart  is  shown  in  Figure  6. 

The  algorithm  is  used  to  monitor  the  wave  while  waiting  for  it  to 
decay  (Figure  4:  prior  to  I  and  between  IV  and  V),  and  to  collect 
the  data  (Figure  4:  VI)  for  subsequent  analysis.  The  routine 
will  provide  the  approximate  location  of  the  extrema  (peaks)  in 
real  time  at  a  scan  rate  of  up  to  75  points  per  second.  If  a  scan 
rate  faster  than  75  points  per  second  is  required,  the  maxima  and 
minima  are  located  after  the  data  have  been  collected  and  before 
initiation  of  the  next  wave. 

In  order  to  digitise  the  signal  efficiently,  the  scan  rate 
(digitisation  rate)  S,  must  be  chosen  to  match  the  characteristics 
of  the  oscillations.  The  optimum  scan  rate  is  a  function  of  the 
period  (P)  of  the  oscillation,  the  number  of  data  points  (N)  col¬ 
lected  per  wave,  the  time  required  for  the  oscillations  to  decay 
to  a  specified  limit,  and  the  method  of  analysis  used.  The  scan 
rate  corresponding  to  40  points  per  cycle  (S  ■  40/P)  provides  an 
adequate  representation  of  the  oscillations  for  most  data  reduc¬ 
tion  methods  (see  later) .  A  rough  estimate  of  the  period  is  ob¬ 
tained  from  the  first  quarter  cycle  after  initiation,  and  the  scan 
rate  is  adjusted  accordingly.  To  locate  the  peaks,  an  interval 
consisting  of  1  +  4R  (where  R  is  a  function  of  scan  rate,  usually 
equal  to  4)  data  points  moves  along  as  the  data  are  acquired,  and 
the  local  maxima  and  minima  are  located  by  determining  whether  the 
center  datum  point  of  the  interval  is  greater  than  (for  a  maximum) 
or  less  than  (for  a  minimum)  both  the  first  and  last  data  points 
of  that  Interval.  As  soon  as  this  set  of  criteria  is  met  the 
center  datum  point  is  stored:  the  next  peak  is  then  sought.  The 
reason  for  using  more  than  three  consecutive  data  points  is  to 
insure  that  a  noisy  signal  does  not  simulate  a  maximum  or  minimum. 
This  method,  although  quite  crude,  is  much  faster  than  one  which 
Involves  taking  a  derivative  of  the  data  to  locate  the  peaks. 

After  all  the  data  points  have  been  collected,  an  approximate 
determination  of  the  peaks  is  made  by  searching  for  tlv:  maximum  or 
minimum  among  the  data  points  within  each  of  the  intervals  in 
which  a  maximum  or  minimum  was  detected. 

Data  Reduction.  The  oscillatory  motion  of  a  freely  moving 
torsion  pendulum  has  been  described  by  an  equation  of  motion  (4): 

I  da8/dt2  +  ndyn  dO/dt  +  Gdyn6  -  0  (5) 

where  I  is  the  moment  of  inertia,  is  the  dynamic  viscosity, 
Gdyn  is  the  elastic  shear  modulus.  6  is  the  angular  deformation, 
ana  t  is  the  time.  The  solution  is  a  damped  sine  wave: 

6  ■  6oexp(-at)cos(ut  +  4)  (6) 


where  0O  is  s  constant;  a  is  the  damping  coefficient* 


“  ‘  "iyn/JI  *  <7) 

u>  is  the  natural  angular  frequency  (radians /sec)* 

.  -  (£f)  -  :  («) 

and  $  is  a  phase  angle.  The  shear  aodulus*  G'»  and  lose  modulus, 
G",  can  be  derived  from  information  in  the  wave: 

G'  -  Kl(ti)2  +  a2)  (9) 

and  G"  -  2Klctu  (10) 

where  K  is  a  geometric  constant. 


Peak  Finding  Method.  Since  the  approximate  location  of  the 
peaks  has  already  been  determined*  the  data  points  about  each  peak 
are  fitted  to  a  quadratic  equation 

6  -  a  +  bt  +  ct2  (11) 

by  a  least  squares  method.  The  optimum  number  of  data  points  to 
be  used  in  fitting  the  quadratic  equation  to  the  data  has  been 
determined  to  be  those  in  the  interval  ♦  0.2*  (7).  Since  the  data 
were  obtained  at  a  scan  rate  such  that  40  points  per  cycle  were 
collected,  the  number  of  points  used  for  fitting  a  quadratic  is 
(0.4*/2*)40  ■  8;  because  the  calculations  require  an  odd  number  of 
data  points,  9  data  points  are  used.  Linear  least  squares  fitting 
of  the  experimental  data  points  to  the  quadratic  equation  requires 
minimisation  of  the  summation  of  residuals  * 


Q  -  I  (f*  -  8.)2  (12) 

i-1  1  1 

where  f^  -  a  +  bt^  +  ct^2  (13) 

and  et  are  experimentally  observed  data  at  tines  t,.  From  the  re¬ 
quirement  that  Q  must  be  minimised* 

U  -  2  l  (f4  -  e.)  (rr)  -  0  (k  -  1  to  3)  (14) 

9  x  i-1  1  1 

where  Aj  -  a,  Aj  -  b,  and  A3  ■  c; 

this  results  in  a  sat  of  three  linear  equations  written  in  matrix 
form: 


(15) 


If  the  time-axis  data  are  offset  ao  that  the  central  datum  point 
la  zero,  the  odd  powered  summations  are  identically  equal  to  zero, 
thus  simplifying  the  matrix.  The  solution  to  this  set  of  equa¬ 
tions  provides  the  parameters  of  the  quadratic  equation.  The  best 
estimate  of  the  peak  position  is  obtained  from  the  first  deriva¬ 
tive 


||  «  b  +  2ct  «  0,  (16) 

U7) 

K  l  4  t.2 

ep"*  +  b<-£)  ♦«(-&)  *  (18) 

This  procedure  is  performed  for  the  first  minimum  and  the  follow¬ 
ing  maximum,  as  well  as  for  the  last  pair  (the  selection  of  which 
depends  on  the  damping),  and  the  period  is  calculated  by  dividing 
the  elapsed  time  between  the  maxima  by  the  number  of  cycles.  The 
logarithmic  decrement  is  obtained  from  the  relation 


6  — 0 

6  -  (A)lntfe*-gg— )] 
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(i  -  3,  5,  7  ...)  (19) 


where  6^  is  the  amplitude  of  the  ith  extremum. 

Least  Squares  Method  (8).  A  torsion  pendulum  specimen  has  a 
tendency  to  change  its  rotational  orlentetion  during  the  course  of 
an  experiment  due  to  an  uneven  distribution  of  stresses  caused  by 
volume  expansion  end  contraction.  This  results  In  a  drift  in  the 
baseline  of  the  wave  signal  which  can  be  represented  by 

8  ■  0Qexp(-at)cos(u)t  +  ♦)  4  Bt  +  C  (20) 

where  B  is  the  drift  coefficient  and  C  is  the  offset.  The  corres¬ 
ponding  differential  equation  can  be  written  as 
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+20^+  (a2  +  u>2)8  -  C(a2  +  w2)  -  2aB  -B(a*+w2)t  ■  0 

(3 

which  may  be  simplified  to 

D"0  +Aljf  +  V+A3t+A4-°-  <J 

Ak  (k  -  1  to  4)  ere  the  parameters  fitted  by  a  linear  lc  ,st 
squares  analysis  to  determine 


u>  -  (A2  -  2 1**  -  ~  .  (24) 

The  derivative  values  of  8  at  any  point  i  are  calculated  numeri¬ 
cally  from  a  quadratic  equation  which  uses  five  consecutive  points 
to  obtain  the  first  and  second  derivatives: 

de.  -26.  ,  -  8.  .  +  6. +  28.  . 

-Jl  -  r _ izZ _ _ £±1 _ ktl}  05) 

Hi-  l  Iflh  J 


d2«i  ( 

dP"  '  t 


20i-2  “  0i-l  “  29i  "  6i+l  +  20i+2-\ 


where  h  is  the  time  interval  between  data  points. 

The  linear  least  squares  fitting  of  n  experimental  data 
points  to  the  differential  form  of  the  equation  of  motion  Involves 
minimization  of  the  summation 

Q  -  l  (f.  -  D.)2  (27) 

i-1  1 

where 

d2e  de. 

fi  "  dtj"  +  A1  dt^*  A20i  +  A3ti  +  A4  <28) 

is  calculated  from  experimental  data  and  D,  is  identically  zero  by 
definition.  From  the  requirement  that  Q  is  minimised* 


S.  • 2  A  Vir  >  ■  0 


(k  •  1  to  4). 


This  set  of  linear  equations  can  be  written  in  matrix  notation: 
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When  this  expression  is  solved  for  A^,  the  values  of  a  and  u  are 
obtained  (Equations  23  and  24). 

Non-Linear  Least  Squares  Method  (9).  Assuming  that  (from 
equation  20) 

6  ■  Ajexp (-A^t) cos (A^t  +  A^)  +  A^t  +  A^  (31) 

where  now  A.  -  6  ,  A2  *  a,  A_  -  2tt/P,  A4  •  $,  A5  ■  B  and  A,  *  C 
is  an  adequate  representation  of  the  solution  to  the  equation  of 
motion  of  a  torsion  pendulum,  the  parameters  A^  (k  *  1  to  6)  can 
be  determined  by  fitting  the  data  (6*,  ti;  1  *  1  to  n)  to  the 
solution. 

If  the  values  of  the  parameters  Afc  were  known,  it  would  be 
possible  to  evaluate 

-  Ajexp(-  A2ti)cos(A3ti  +  A^)  +  A5t±  +  Afi  (32) 

for  each  t^  to  obtain  a  set  of  "true"  residuals 

r±  -  -  eA  (i  -  1  to  n).  (33) 

A  "true"  residual  would  represent  the  difference  between  the 
actual  function  value  at  t*  and  the  empirical  value  6*.  These 
"true"  residuals  cannot  be  calculated  because  the  actual  values  of 
the  parameters  Afe  are  not  known. 

However,  initial  estimates  of  the  parameters  A£  can  be  ob¬ 
tained  from  other  methods,  or  a  previous  wave,  and  "computed" 
residuals  can  be  calculated: 

R^ ■  A^exp(-A*tj)cos(Ajt^  +  A^)  +  A^t^  +  “  ®i  (34) 

(i  »  1  to  n) 
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Improved  estimates  of  the  parameters  A^  can  be  obtained  by  a  dif¬ 
ferential  correction  technique  based  on  least  squares,  provided 
that  the  estimates  a£  are  sufficiently  close  to  the  actual  values 
of  the  parameters  Aj<  to  lead  to  convergence  of  the  method.  This 
differential  correction  technique  can  be  derived  by  first  expand¬ 
ing  the  function  about  A&  using  a  linear  Taylor  series  expansion 
of  the  form 


f(t^»  Aji  A2*  •••»  A^)  *  f (t^i  A^i  A^»  ...»  Aj) 


+  3Ax(A1  “ 


Al}  +  3 A 2<A2  " 


A2> 


+  11  (A  - 
3A6  6 


a;> 


(35) 


so  that  a  relation  between  the  r^  and  Ri  can  be  obtained.  This 
relation  can  be  found  by  evaluating  the  equation  at  each  value  of 
t£  and  subtracting  from  both  sides  of  the  equation.  Using  the 
definitions 


6Ak  *  \  “  \ 


and 


3f 


3f , 


- i-  .  -L 

3Ak  3\ 


|t  *  V  Ak  «  Ak; 
the  result  can  be  written  in  the  form 


(36) 


(37) 


o  o  o 

^ •••*  A^)  *  6^  ■  f (t^»  •  ••*  A£) 


'6'  wi 

+  (3^)  «*!  *  +  <^> 


3f . 


2*  •  •  • »  "g' 

(i  -  1  to  n)  (38) 


The  desired  relation  between  the  r^  and  can  then  be  found  by 
substituting  the  expressions  for  r^  and  R*: 


3fi 

Ri  +  (1^)SA1  + 


3fi 

+  <3a7>6V 

O 


(i  -  1  to  n).  (39) 


This  relation  can  be  used  to  compute,  from  a£,  a  set  of  parameters 
Afe  that  minimizes  the  sum  of  the  squares  of  the  "true"  residuals 
rj,  i.e.« 


n  y  n  3f . 

Q  -  I  *7  -  I  1*,  ♦  Car)6Ai  + 

1*1  *  i*l  *  * 


3fi  2 

**  +  (3A6)6A61  ’ 


(AO) 


The  function  Q  has  a  minimum  value  when  all  of  its  partials  with 
respect  to  the  6Ajc  are  simultaneously  zero: 
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(45) 


(46) 


Ifl 

BA, 


5 

3A, 


BA. 


BA, 


-A1tiexp(-A2ti)«in(A3t1  +  A4), 
-A1exp(-A2t1)sin(A3t1  +  A4), 


(47) 

(48) 

(49) 


The  solution,  fA^,  of  this  set  of  normal  equations  is  a  first- 
order  approximation  of  the  changes  in  A^  required  to  obtain  the 
parameters  Afc.  If  any  |6Afc.|  >  c  (error  limit),  Afc  is  replaced  by 
Afc  -l-  5Ak  and  the  entire  differential-correction  procedure  is  re¬ 
peated  using  these  new  estimates. 


Fourier  Transform  Method.  Another  method  of  data  reduction 
is  to  take  a  fast  Fourier  transform  (FFT)  of  the  wave  (10) .  As  indi¬ 
cated  in  Figure  7,  the  Fourier  transform  of  a  damped  sine  wave 
with  a  single  frequency  is  a  single  maximum  in  the  frequency  . 
domain  at  the  frequency  of  the  oscillation.  The  amplitude  (H)  of 
the  transformed  data  as  a  function  of  angular  frequency  (u)  is 
given  by  (11) 


6  [a2 (a2+w24w2)2  +  w2(o2+w2-u2)2]1* 

(a2+ui2-u2)2  +  (2au>)2 
o 


(50) 


0  is  the  initial  amplitude,  a  is  the  damping  coefficient. 


where  8 

and  w0  is  the  natural  angular  frequency  of  the  oscillation, 
amplitude  of  the  peak  is  given  by 


max 


8 


a2  + 


(( 


or  +  4u>’ 


§■>]*. 


e 

0 

2a 


(a  «  uQ) 


The 


(51) 


from  which  a,  the  damping  coefficient,  is  obtained. 

Discussion.  The  four  methods  of  data  reduction  were  used 
to  analyze  the  raw  data  of  the  same  TBA  specimen  during  a  slow 
(0.25°C/mln)  temperature  scan  (Figure  8).  A  comparison  of 
the  spectra  indicates  that  they  all  gave  similar  results  over 
the  range  of  period  (0.3  to  1.8  sec)  and  logarithmic  decrement 
(0.01  to  1.08)  encountered  in  the  experiment.  (The  automated 
torsion  pendulum  has  been  used  to  reduce  data  with  a  range  of 
0.1  to  15  sec.  for  the  period,  and  0.001  to  4.0  for  the 


logarithmic  decrement).  In  Figure  8  there  is  no  appreciable  dif¬ 
ference  in  the  relative  rigidity ,  but  in  the  logarithmic  decrement 
the  non-linear  least  squares  reduction  method  produces  the  smooth¬ 
est  results,  followed  by  the  peak-finding  method.  The  linear 
least  squares  and  Fourier  transform  methods  have  considerable 
scatter.  The  linear  least  squares  method  also  results  in  a  small 
systematic  difference  in  the  logarithmic  decrement  from  the  other 
three . 

Although  the  peak-finding  method  is  the  simplest  way  of  de¬ 
riving  the  period  and  logarithmic  decrement  from  the  raw  data,  it 
has  some  limitations  in  that  at  least  2.25  cycles  of  oscillation 
are  required  in  order  to  do  the  calculation.  This  is  a  problem 
when  the  system  approaches  critical  damping  conditions.  It  also 
is  difficult  to  calculate  the  logarithmic  decrement  at  very  low 
damping  (when  the  peak  amplitude  changes  only  slightly  during  the 
time  data  is  collected)  due  to  the  resolution  of  the  digitizing 
voltmeter.  Of  the  four  methods  discussed,  the  peak  finding  method 
is  the  most  sensitive  to  the  scan  rate,  since  the  number  of  data 
points  about  each  peak  that  are  fitted  to  the  three-parameter 
quadratic  equation  should  be  at  least  nine  as  discussed  earlier. 
Therefore  the  scan  rate,  which  depends  on  the  estimate  of  the 
period,  needs  to  be  quite  close  to  40/P,  where  P  is  the  actual 
period.  Also,  a  systematic  error  is  Introduced  when  this  method 
is  used,  because  the  position  of  the  peaks  is  a  function  of  the 
damping  as  well  as  the  period.  As  can  be  seen  in  Figure  9,  the 
peaks  shift  to  shorter  times  as  the  damping  increases. 

The  least  squares  fitting  method  overcomes  limitations  of  the 
peak  finding  method,  but  introduces  some  errors  of  its  own.  It  is 
not  as  sensitive  to  the  scan  rate,  as  long  as  S  »  2/P  (the 
Nyquist  frequency) ,  and  its  accuracy  increases  with  increasing 
number  of  data  points.  It  was  found  empirically  (8)  that  the 
error  was  reduced  if  an  integral  number  of  cycles  was  used  in  the 
analysis.  Some  error  is  Introduced  due  to  the  fact  that  first  and 
second  derivatives  of  the  raw  data  have  to  be  taken. 

The  use  of  the  non-linear  least  squares  method  does  not  re¬ 
quire  any  derivatives,  but  needs  an  initial  estimation  and  takes 
more  time  to  compute,  since  several  Iterations  (usually  3  or  4) 
are  necessary  to  reduce  the  difference  between  the  estimated  and 
calculated  values  of  the  damping  coefficient  to  within  0.1Z.  But 
since  this  method  only  requires  between  100  and  150  data  points 
without  a  loss  in  accuracy  compared  to  as  many  as  1000  for  the 
peak-finding  and  least  squares  methods,  the  scan  rate  can  be  re¬ 
duced  as  much  as  90Z  and  the  time  required  for  the  calculations  is 
reduced  to  the  order  of  a  minute. 

The  Fourier  transform  method  requires  a  minimum  of  1024  data 
points  to  provide  enough  resolution  to  calculate  the  damping  co¬ 
efficient.  The  FFT  of  1024  data  points  takes  approximately  a  min¬ 
ute  with  the  HP  9825B  computer,  so  this  constitutes  a  practical 
limit  in  resolution  due  to  computer  memory  size  and  time  consider¬ 
ations.  So  as  not  to  Introduce  error,  the  damped  oscillations 


oust  not  be  truncated;  therefore  it  is  laportant  to  adjust  the 
scan  rate  and  the  number  of  data  points  so  that  the  entire  wave  is 
collected.  The  major  problem  with  using  the  FFT  method  la  the 
difficulty  in  obtaining  accurate  values  of  a;  the  curve  is 
Laurentian,  and  hence  its  amplitude  at  the  maximum  is  difficult  to 
obtain.  One  way  around  this  is  to  use  a  curve  fitting  procedure, 
but  then  there  is  no  advantage  in  using  this  method.  Some  alter¬ 
natives  may  be  to  use  a  larger  and  faster  computer,  or  a  dedicated 
microprocessor  such  as  the  spectrum  analyzer  (HP  3582A) ,  which  can 
compute  the  FFT  in  real  time.  A  practical  feature  of  the  FFT  is 
in  the  display  of  the  transformed  data;  any  non-homogeneity  of  the 
signal  due  to  other  modes  of  motion  will  appear  as  secondary  peaks, 
and  so  this  method  serves  as  an  excellent  way  to  monitor  the 
oscillations. 

Although  comparison  of  the  four  methods  shows  that  the 
smoothest  reduced  data  for  the  given  experiment  were  obtained 
using  the  non-linear  least  squares  method,  the  ultimate  quality 
depends  on  the  quality  of  the  sensor  signals  of  the  experiment. 
Published  superior  TBA  spectra  obtained  using  the  linear  least 
squares  method  (8) ,  and  the  peak-finding  method  using  an  analog 
computer  (2),  were  presumably  the  consequence  of  a  better  basic 
experiment~~than  the  one  used  in  this  report  to  compare  (as  in 
Figure  8)  the  data  reduction  methods. 

Calibration 


A  calibration  wire  whose  shear  modulus  is  known  can  be  used 
to  determine  the  moment  of  inertia  of  the  pendulum  assembly,  so 
that  quantitative  meaaureDents  of  the  dynamic  mechanical  proper¬ 
ties  of  specimens  can  be  made.  The  shear  modulus  of  the  calibra¬ 
tion  wire  la  obtained  by  msaaurlng  the  period  of  oscillation  of  a 
simple  torsion  pendulum  consisting  of  an  aluminum  rod  suspended  by 
the  wire.  The  moment  of  inertia  of  this  system  is  given  by 


I 


(52) 


where  m  is  the  mass,  r  is  the  radius  and  l  la  the  length  of  the 
rod. 

The  she*?  modulus  G’  of  a  wire  is  given  by 


C' 


(53) 


dmr«  L  is  the  length  and  R  is  the  radius  of  the  wire.  Vlth  the 
calibration  wire  (^oee  shear  modulus  was  determined  to  be  9.789  * 
10 11  dyne /cm*)  as  a  specimen  in  the  automated  torsion  pendulum, 
the  moment  of  inertia  of  the  pendulum  that  was  used  routinely  in 
the  subsequent  experiments  was  determined  to  be  138.7  g-ce2.  The 
shear  modulus  of  a  film  of  known  dimensions  can  then  be  calculated 


•/ 
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from  the  period  end  logarithmic  decrement  using  the  equation  (4): 


4*aIL 


Of*-  (1  +  '  ijs 


where  K  Is  a  form  factor: 

ab* 


H  -  -=|-  (1  -  0.63  b/a). 


(54) 


(55) 


a  is  the  width,  b  is  the  thickness  (b  <  a/3),  L  is  the  length,  m 
is  the  mass  supported  by  the  specimen  and  g  is  the  gravitational 
constant. 


Comparison  of  Torsion  Pendulum  and  TBA 

A  film  of  an  amine-cured  epoxy,  Epon  828  (Shell) /PACM-20 
(DuPont)  with  Tg*>  *»  166°C,  was  cured  by  heating  it  to  250°C  in  a 
helium  atmosphere.  The  dynamic  mechanical  spectrum  of  this  film 
is  shown  in  Figure  5.  For  comparison,  the  corresponding  spectrum 
of  a  specimen  consisting  of  a  multifllamented  glass  braid  impreg¬ 
nated  with  the  uncured  resin  and  cured  in  the  TBA  apparatus  by 
heating  it  to  200SC  under  helium  atmosphere  is  shown  in  Figure  10. 

It  has  been  reported  (12)  that  the  shear  modulus  as  measured 
by  the  torsion  pendulum  and  TBA  should  differ  only  by  a  multipli¬ 
cative  constant  below  Tg  and  the  logarithmic  decrement  should  be 
identical.  Although  the  spectra  of  Epon  828/PACM-20  obtained  by 
torsion  pendulum  and  torsional  braid  analysis  show  transitions  at 
the  same  temperature  (glass  transition  at  166°C  and  a  secondary 
sub-glass  transition  at  -28°C),  the  results  indicate  that  the 
actual  modulus  and  logarithmic  decrement  cannot  be  compared  quanti¬ 
tatively.  In  Figure  11  the  relative  rigidity  (TBA)  has  been 
shifted  vertically  for  comparison  with  the  torsion  pendulum  data; 
a  vertical  shift  on  a  logarithmic  scale  is  equivalent  to  multi¬ 
plying  by  a  constant.  It  is  evident  by  comparing  the  curves  in 
Figure  11  that  there  is  only  a  qualitative  correlation  between  them. 

Conclusions 


The  automation  of  the  torsion  pendulum  utilising  a  desktop 
computer  eliminates  the  tedious  data  analysis  previously  associated 
with  that  technique.  Any  one  of  four  data  reduction  methods  can 
be  used;  the  experimental  conditions  will  determine  which  is  the 
optimum  one  to  employ.  The  torsion  pendulum  technique  provides 
quantitative  values  of  shear  modulus  and  logarithmic  decrement  and 
in  the  torsion  braid  mode  provides  a  qualitative  analysis  of 
■aterials,  especially  in  the  llquid-to-solld  transition  region. 

In  addition  to  providing  the  capability  of  using  any  one  of  four 
data  reduction  techniques,  the  computer  has  the  advantage  of  stor¬ 
ing  the  data  on  magnetic  tape,  where  it  is  available  to  be 
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accessed  for  further  computation  or  to  be  plotted  in  whatever  mode 
Is  most  suitable.  Since  the  computer  is  easily  programmable,  the 
software  can  readily  be  adapted  to  consider  other  variables  or  to 
control  the  experiment  in  other  ways. 
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Figure  Captions 

Figure  1.  Automated  torsion  pendulum:  ecbeaatic.  An  analog 

electrical  signal  is  obtained  from  passing  a  light  beam 
through  a  pair  of  polar izera,  one  of  which  oscillates 
with  the  pendulum.  The  pendulum  is  aligned  for  linear 
response  and  initiated  by  a  computer  that  also 
processes  the  damped  waves  to  provide  the  elastic 
modulus  and  mechanical  damping  data  which  are  plotted 
on  an  XYY  plotter  versus  temperature  or  time. 

Figure  2.  Automated  torsion  pendulum:  system  schematic  for 
Interfacing  with  a  digital  computer.  The  torsion 
pendulum  has  been  interfaced  with  a  digital  desktop 
computer  (HP-9825B) .  The  motors  which  align  the 
specimen  and  initiate  the  waves  are  under  computer 
control.  The  wave  and  amplified  analog  thermocouple 
signals  reach  the  computer  digitized  via  a  digital 
voltmeter  (HP-3437A) .  The  scanner  (HP-3495A)  super¬ 
vises  the  I/O  activity.  Upon  receiving  the  digitized 
raw  data  the  computer  calculates  the  frequency  and 
damping  parameters,  and  plots  the  dynamic  mechanical 
properties  of  the  specimen  as  a  function  of  tempera¬ 
ture  and  time. 

Figure  3.  Automated  Torsion  Pendulum:  the  pendulum  is  housed  in 
the  cabinet  at  the  left;  the  oven  is  separated  from 
the  optical  transducer  by  an  insulated  3/4  inch 
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horizontal  aluminum  plate.  The  temperature  controller* 
digital  voltmeter, scanner,  and  computer  are  in  the 
rack  at  the  right.  The  atmosphere  control  panel  and 
liquid  nitrogen  container  are  shown  in  the  background. 
One  of  the  authors  (JBE)  la  seated  at  the  console. 
Figure  4.  Automated  torsion  pendulum:  control  sequence. 

I)  Previous  wave  decays,  drift  detected  and  correc¬ 
tion  begins.  Zl)  Reference  level  of  polarizer  pair 
reached.  Ill)  Wave  initiating  sequence  begins. 

IV)  Decay  of  transients.  V)  Free  oscillations  begin. 
VI)  Data  collected.  VII)  Control  sequence  repeated. 
Figure  5.  Dynamic  mechanical  spectrum  (torsion  pendulum)  of  a 
cured  film  of  Epon  828/PACM-20.  Both  the  shear 

modulus  G'  (O)  and  its  approximation  G’  (0)  are 

% 

plotted  on  the  upper  curve;  the  lower  three  curves 
are  loss  modulus  G"  (o) ,  logarithmic  decrement  &  (*) , 
and  damping  coefficient  a  (+). 

Figure  6.  Flow  diagram  of  data  collectlon/peak-flnding 
algorithm. 

Figure  7.  Fourier  Transform  Method.  The  Fourier  transform  of 
an  exponentially  damped  sine  wave  of  period  P  and 
damping  coefficient  a  is  a  single  maximum  at  the 
oscillation  frequency  whose  amplitude  is  Inversely 
proportional  to  the  damping  coefficient. 
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Figure 


Figure 


Figure 


Figure 


8.  Dynamic  mechanical  spectrum  (TBA)  of  Epon  828/PACM-20 
in  which  the  relative  rigidity  and  logarithmic 
decrement  have  been  calculated  by  four  methods: 
peak-finding  [PKF  (□)],  linear  least  squares 

[LSQ  (*)  .  non-linear  least  squares  [NLSQ  (+)],  and 
fast  Fourier  transform  [FFT  (o)].  For  clarity  the 
LSQ,  NLSQ  and  FFT  data  have  been  displaced  vertically 
in  equal  increments  from  the  PKF  data. 

9.  Damped  Sine  Wave.  Error  in  measuring  period  by  peak- 

finding  method:  the  peaks  of  an  exponentially 
damped  sine  wave  of  single  frequency  (0.5  Hz)  shift 
to  shorter  times  with  Increasing  damping  coefficient 
(a  values:  —  0.05,  —  0.5,  -.-  1.0,  - - 1.5). 

10.  Dynamic  mechanical  spectrum  (TBA)  of  a  cured  composite 
specimen  (glass  braid  impregnated  with  Epon  828/PACM- 
20  resin) . 

11.  A  comparison  of  torsion  pendulum  data  obtained  using 
a  film  [G'  (Q),  A  (*)]  and  TBA  data  obtained  using 

a  supported  specimen  [relative  rigidity  (0),  A  (o)]. 
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